Superconducting Vortex Lattices for Ultracold Atoms 
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The ability to trap and manipulate ultracold atoms in lattice structures has lead to a remarkable 
experimental progress to build quantum simulators for Hubbard models. A prominent example 
is atoms in optical lattices where lasers are used to create lattices with spacing set by the laser 
wavelength as well as to control and measure the many-body states. In contrast, here we propose and 
analyze a nanoengineered vortex array in a thin-film type-II superconductor as a magnetic lattice for 
ultracold atoms. This proposal addresses several of the key questions in the development of atomic 
quantum simulators. By trapping atoms close to the surface, tools of nanofabrication and structuring 
of lattices on the scale of few tens of nanometers become available with a corresponding benefit in 
energy scales and temperature requirements. This can be combined with the possibility of magnetic 
single site addressing and manipulation together with a favorable scaling of superconducting surface 
induced decoherence. 



Optical lattices allow to confine ultracold neutral 
atoms in highly controllable periodic potentials de- 
signed by the interference pattern of counterpropagat- 
ing lasers [1, 2]. This possibility has triggered sev- 
eral research fields during the last decade, in particular, 
quantum simulation of many-body physics [3, 4], quan- 
tum computation [5], and quantum metrology [6]. Re- 
markably, the quantum simulation of Bose- and Fermi- 
Hubbard hamiltonians holds the promise to be the first 
technology that can solve problems not accessible to clas- 
sical computers [7]. An important challenge to achieve 
this goal is the development of better cooling schemes in 
order to reduce the entropy of the simulator [3]. Effec- 
tively, this could be achieved by boosting the energy scale 
of the physical parameters of the quantum simulator by 
scaling down the lattice potential, namely by reducing 
the distance between lattice sites. Due to the diffraction 
limit, this requires to trap atoms near a surface without 
adding new sources of decoherence, which results in a 
challenging trade-off, see for instance the recent proposal 
based on nanoplasmonic optical lattices [8]. However, 
placing atoms close to surfaces also makes available the 
tools of nanofabrication to engineer arbitrary (periodic 
and non-periodic) atomic lattice structures. 

Here, we propose and analyze a new approach to trap 
and manipulate ultracold neutral atoms in arbitrary lat- 
tice potentials based on using a magnetic nanolattice gen- 
erated by a controlled array of superconducting vortices 
in thin films. While magnetic traps for ultracold atoms 
created by macroscopic currents flowing in type-H su- 
perconducting materials (the ones that contain vortices) 
have been theoretically proposed and experimentally re- 
alized [9-14], here, we aim at something radically dif- 
ferent and more challenging, namely to trap single atoms 
using the field created by few controlled superconducting 
vortices. That is, we propose to exploit the nanoscience 



and technology developed to manipulate superconduct- 
ing vortices in thin films [15] for atomic physics. This 
technology has the goal of controlling the electronic prop- 
erties of superconducting devices that are strongly af- 
fected by the presence and motion of vortices. In partic- 
ular, today superconducting vortices can be positioned 
more or less at wish, even in complex structures, by ar- 
tificially pinning them in nanoengineered arrays of, for 
instance, completely etched holes (antidots) of various 
sizes and shapes, see [15] and references therein. Our 
proposal hints at the possibility to fabricate and struc- 
ture magnetic lattices at the fundamental length scales 
associated to superconducting vortices, which can be of 
few tens of nanometers, and hence render magnetic lat- 
tices to an unprecedented hitherto parameter regime [16- 
19]. Moreover, the combination of all-magnetic trapping 
and manipulation and superconducting surfaces leads, in 
principle, to very favorable scalings on surface-induced 
decoherence, as discussed below. 

Superconducting Vortex Lattice 

We consider a superconducting film of thickness d < X, 
where A is the London penetration depth. The film is a 
type-H superconductor, i.e., A/^ > 1/a/2, where ^ is the 
coherence length whose value is typically between few 
to tens of nanometers [20]. In thin films, the effective 
penetration depth is given byA = A^/(i>A>(i>^, 
which hence can potentially be as small as few tens of 
nanometers. The upper side of the film is situated at 
the X — y plane with z = and contains a periodic 
array of superconducting vortices distributed in a Bra- 
vais lattice R = niai + n2a2, where ai(2) are the lattice 
primitive vectors and ni(2) range through all integer val- 
ues. The density of vortices is 1/a^ where = |ai x a2| 
is the area of the primitive cell. Each vortex creates a 
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FIG. 1: Superconducting vortex nanolattice. a) An array of completely-etched holes (antidots) with a commensurate 
number of vortices creates a magnetic field that can be used to trap neutral atoms in a two-dimensional magnetic lattice with 
the same lattice geometry as the vortex lattice, b) The distance between vortices is given by a and atoms are trapped at 
a distance zq from the upper side of the film. The magnetic field generated by a vortex can be approximated by the field 
generated by a monopole situated at z = —dm A of charge 2$o- The trapping along the z-axis is achieved by applying a 
perpendicular and homogeneous magnetic field of strength Bi . A time-dependent bias field parallel to the film Bm (t) is applied 
to create a time-averaged orbiting potential which is employed to circumvent Majorana losses. 



magnetic flux given by <l>o = /i/(2e), where h = 27rh is 
Planck's constant and e the charge of the electron. In su- 
perconducting films with a low intrinsic pinning, a trian- 
gular lattice minimizes the vortex- vortex repulsive inter- 
action [20]. In this case, the maximum vortex density is 
limited to 1/^^, which corresponds to the situation when 
the cores of the vortices overlap. Other lattice geometries 
can be considered by nanoengineering a periodic array of 
artificial pinning centers consisting of antidots of radius 
R > ^ [15]. The artificial pinning sites attract the vor- 
tices very strongly [21, 22] and, as a result, the underlying 
symmetry of the pinning array dictates the symmetry of 
stable vortex patterns. The interplay between the inter- 
action of vortices with the pinning lattice and the vortex- 
vortex repulsion leads to interesting phase diagrams for 
the crystallization of vortices [27]. In the following, we 
assume that the vortices stabilize at the pinning lattice 
if the density of antidots does not exceed 1/a^ < 
see Supplemental Information (SI) for further discussion. 
The vortex lattice can be prepared by cooling the film in 
the presence of an homogeneous perpendicular matching 
field of strength Bg, = ^q/o? such that a commensurate 
number of superconducting vortices is created in the anti- 
dot lattice. At the final temperature below (which can 
be of a high-Tc material) the applied field is switched-off 
and the vortices remain. 



Magnetic field generated by the vortex lattice 

Our proposal consists in using the magnetic field cre- 
ated by the vortex lattice above the film Bv(r, z > 0) 
to trap neutral atoms in a two-dimensional magnetic 



nanolattice with a geometry dictated by the nanoengi- 
neered antidot lattice, see Fig la. Hence, we need to 
calculate Bv(r,2: > 0), which for an infinite arbitrary 
array of superconducting vortices is a complicated task. 
For a triangular lattice, this has been discussed within 
the Ginzburg-Landau theory [23], but its expression is 
cumbersome. In the following, we will make use of a sim- 
plified model that allows us to obtain Bv(r, z > 0) very 
accurately with a simple expression. First, we consider 
the London limit A ^ ^, such that the currents and field 
generated by a vortex in a thin film can be analytically 
obtained [24], see SI. Furthermore, as discussed in [25, 26] 
and in the SI, the field above the film created by a single 
vortex situated at R, which we denote by BK(r, z > 0), 
can be very well approximated by the one created by a 
magnetic monopole situated at 2: = — <im ^ —A and with 
a charge 2^o such that a flux ^0 passes through the semi- 
circle of radius [|r-R|2-h(z-h(im)^]^/^ That is, we will 
use 



BR(r,^>0) 



^0 



(r-R,2; + dm) 



2^ [|r-R|2 + (z + d^)2]^/2' 



(1) 



Using this approximation, see Fig. lb, we can calcu- 
late the field By generated by an arbitrary vortex lat- 
tice {R} by solving the Poisson equation V^(/)(r, z) = 
— p(r, z), where Bv(r, z) = — V(/)(r,z), p(r, z) = 
2<l>o(^ {z + dm) X]r ^ (^^ ~ 1^)5 (5() is the Dirac delta 
function. This can be analytically solved by Fourier 
transformation into the reciprocal space, where in this 
case p(r, z) = 2Bq8 {z ^ dj^)^j^e~'^^'^ . The sum is 
over all reciprocal lattice vectors K and we have defined 
^0 = ^o/o?' As shown in the SI, the solution of the 
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Poisson equation leads to 



B^''>{r,z) = Boe-'^'g,^y){r,z), 



BUr,z) = Bo[l 



(2) 



where we have defined = 2k{z + d,^) > 2kd,^ = Amin, 
k = Tr/a, and, 



5.(,)(r,. > 0) ^ ^ ^sin(K.r)e^^(-l^l/'^), 



K/0 



,(r,z>0)= ^cos(K-r)e^^(i-l^l/^). 



(3) 



The X and y components of By are equal zero on top 
of the vortices, namely at r = R. The z component 
is always positive and tends to an homogeneous field of 
strength Bq at long distance from the surface. The non- 
zero value of the field at long distances is a consequence 
of the infinite extension of the plane. Note that general 
non-periodic structures obtained by nanoengineering can 
also be considered. 
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Magnetic lattice potential for ultracold atoms 

Alkali metal atoms in low fields of strength < 30 
mT, where the Zeeman shift of hyperfine levels is lin- 
ear, experience a potential of the form Mat(r,2;) = 
/i^^ |B(r, 2:)| [28, 29]. The local field interacting with 
the atoms is denoted by B and will be composed by the 
one generated by the vortices By plus additional bias 
fields, see below. The magnetic dipole moment is given 
by jiimp = rriFgFl^B^ where uif is the magnetic quan- 
tum number, is the Lande ^-factor, and the positive 
number /iB is the Bohr magneton. Thus, low-field-seeker 
states qf'^f > can be trapped at the local minima of 
|B(r, z)\ [28, 29]. The form of the potential assumes that 
the projection of the atomic magnetic moment into the 
instantaneous direction of magnetic field is constant, an 
assumption that we take into account below. 

The field generated by the vortices does not have local 
minima of |Bv| since as discussed in the previous section, 
the 2;-component is always positive. For this reason, we 
propose to add a perpendicular bias field of the form 
Bi = 5i (0,0,-1) and define Biat(r,z) = Bv(r,z) + 
Bi. The field Bi can be considered homogeneous even 
close to the surface provided that a thin film of d < A is 
used. We have validated this assumption by numerically 
calculating the field distribution and induced currents 
in a superconducting disk of finite radius and thickness 
using an energy minimization procedure [30] taking into 
account their kinetic and magnetic energy terms [31], see 
SI. The strength of Bi is limited by the fact that it should 
not induce vortices in the central region of the film. This 
leads to the condition Bi < B^ -\- minr5Y(r,0), where 
B"^ ^ <l>o/(47rA^). The local minima of |Biat(r, 2:)| are 



FIG. 2: Lattice control with the perpendicular bias 
field, a) Relation between the applied field Bi and the di- 
mensionless parameter A = 27r(2;o -\-dm)/a. Since zo can only 
be positive, A has a minimum given by Amin = 27rdm/cL. Bi 
is also upper-limited by the field B^ in order not to induce 
vortices. At = there is no confinement along the z- 
axis. In the inset we plot |Biat(x, 2/, ^o)|^ for a square lattice 
and A = 2 for > Bq (deep wells) and Bi < Bq (sharp 
peaks), b) |Biat|^ as a function of x/a at 2/ = and z = zq 
for a square lattice for different A. |Biat|^ is plotted in units 
of Bmax = |Biat(a: = a/2,0,zo)|^- As seen in the plot, the 
smaller the A, the more Fourier components are involved in 
|Biat|^ as a function of x. 



obtained on top of the vortices, r = R, when Bi > Bq, 
at a position zq given by the solution of the equation 
5y(R, zq) = 5i, see Fig. 2a, where we have defined A = 
A^o- In this case, the conditions B^ > Bi > Bq can 
be fulfilled provided A > a. In Fig. 3 we plot |BiatP 
(as the final form the potential depends on the modulus 
square, see Eq. (4)) in the x — z plane (at y = 0). As 
a validation of the model, we plot the same quantity for 
an array of 3 x 3 vortices in a finite plane numerically 
solving the London equation using the method presented 
in [32]. One can observe that there is also a minimum 
of |Biat(r, z)p on top of the central vortex. In the case 
Bi < Bq^ one confines atoms between vortices with a 
much shallower trap, see Fig. 2a. 

It is well known that Majorana losses [33], namely 
spontaneous spin fiips rendering the state of the atom 
into a high-field-seeker, occur when |B| ~ 0. Note that 
the minima of |Biat| correspond to zero-field. In order 
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FIG. 3: Magnetic field distribution, a) A contour plot of 
|Biat|^ at 2/ = is shown for a square lattice as a function of 
xja and = 27r(^ + (im)/a for A = 27r(2;o + c?m)/tt = tt. This 
is done using Eq. (2), which is calculated using the approxi- 
mation of the monopoles and for an infinite lattice. The local 
minima |Biat|^ on top of the vortices are clearly observed, b) 
The same plot as a) is done with a set of 3x3 localized vor- 
tices located on a squared lattice in a thin superconductor. 
|Biat|^ is calculated by numerically solving the London equa- 
tion using the method presented in [32] and using A a. A 
local minimum is also observed on top of the central vortex. 
Note also the boundary effects (the borders of the plane are 
at xja = ±2.45. 



to prevent Major ana losses we suggest to use an effec- 
tive time- averaged, orbiting potential [34] generated by 
adding to Biat an homogeneous time-dependent field par- 
allel to the thin film Bm(^) = ^mIsIucjm^, coscjm^, 0). 
By time- averaging we have that (Bm(^)) = 0, but 
(|BM(t)|) = ^M- Assuming > max |Biat(r, zq)! 
and considering that the total field experienced by the 
atoms is given by B = Biat + Bm, we have (|B|) ^ 
Bm + |BiatP/(25M), which does not contain zero-field 
local minima. Using that the oscillation frequency of the 
atom cjt is much smaller than cjm, and cjm much smaller 
than the Larmor frequency (so that the atomic potential 
depends on the modulus of the magnetic field), namely 
<^ <C cjl = Mtuf^m/^, the effective time- averaged 
magnetic potential for the atoms is given by [34] 

Vlat(r,^)«fi^L + ^|Biat(r,^)p. (4) 
This potential depends on |BiatP, avoids Majorana 



losses, and confines atoms on a magnetic lattice whose 
geometry is dictated by the vortex lattice. Note that the 
bias field B\ can be used control the trapping position in 
the z-axis, see Fig. 2a. This might be used to adiabati- 
cally load the ultracold atoms into the magnetic lattice 
from an external trap. 

Quantum simulation of Hubbard hamiltonians 

A prominent application of the magnetic lattice poten- 
tial Eq. (4) is the possibility to perform quantum simu- 
lation of quantum many-body physics. When loading 
an ultracold gas of neutral atoms into the magnetic lat- 
tice, atoms will be positioned at the local minima of 
the lattice potential. In this situation, atoms can tun- 
nel to neighboring lattice sites and interact on-site due 
to short-ranged collisional interactions [2]. This leads to 
similar rich many-body physics encountered in typical 
condensed-matter physics, where electrons can be mod- 
eled as moving on a lattice generated by the periodic 
array of atom cores. The advantage of quantum simu- 
lators is the high-degree of controllability of purely con- 
servative and defect-free lattice potentials. The hamil- 
tonian describing the interaction of ultracold atoms in 
a lattice is given by a Hubbard model, which is one 
of the most paradigmatic models in condensed-matter 
physics [4]. For bosonic atoms, the Hubbard hamilto- 
nian reads 

H = —t ^ + o^-cij^ + U ^ fiifii^ (5) 

where ai{a\) annihilates (creates) a bosonic localized 
atom on site i, and fii = a\ai. The first term describes 
tunneling between neighboring sites with tunnel coupling 
t, and the second on-site interactions with strength U . 
The particular interest of studying these type of models, 
in particular for spin=l/2 Fermionic atoms [35, 36], lies 
on the fact that in the strong coupling regime U/t ^ 1, 
superexchange processes provide the basic mechanism for 
an antiferromagnetic coupling between spins on neigh- 
bouring sites, which is closely related to studies of high- 
Tc superconductivity within the Hubbard model [37]. 

The strength of the Hubbard parameters and the de- 
pendence with the physical parameters of the supercon- 
ducting vortex lattice (SVL) proposed here can be readily 
obtained in analogy to optical lattices (OLs) [1]. In par- 
ticular, let us consider a square lattice with Bi > Bq and 
A sufficiently large, such that the potential Eq. (4) can 
be approximated to Mat(r, ^) ~ Vb[sin^(A:x) +sin^(%)], 
where Vq = SBq exp[— 2A]/i^^/5M. This potential has 
the same form as the typical one obtained in OLs. The 
distance between the trapped atoms in the SVL is given 
by a, therefore 2a plays the role of the optical wavelength 
in OLs. The role of the laser intensity in OLs, that can be 
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used to modulate the trap depth, is taken in the SVL by 
the strength of the bias field Bi , recall that A depends on 
Bi, see Fig. 2. Therefore, by defining = h'^k'^ /{2ma), 
one can use the same formulas used in OLs to obtain 
the parameters in Hubbard hamiltonians [2, 35]: the 
trap frequency uJt ~ ^/4E^yo/h^ the tunneling rate 
t ^ 2En{Vo/Enf^exp[-2{Vo/Eny/^]/^, and the 

3/4 



on-site repulsion U ^ y S/TrE^Rag/c (Fq/^r) , where 
is the 5-wave scattering length. In Fig. 4, we plot the 
tunneling rate as a function of a/ A for Rb and Li (for 
A = 100 nm and A = 20 nm) and compare it with the 
decoherence rates that we discuss below. We also pro- 
vide typical values for the trap frequencies in the caption 
of the figure. Since half the optical wavelength in some 
OLs can be of 250 nm [1], the case A = 20 nm supposes 
two orders of magnitude larger values of £^r, and there- 
fore, significantly lower temperature requirements for the 
simulation of quantum magnetism in Hubbard hamilto- 
nians [36]. 

The potential Eq. (4) also permits to design dense lat- 
tices with higher Fourier components, something that is 
very challenging to achieve in OLs. Here, by reducing A, 
more reciprocal vectors of different frequencies enter into 
play, see Fig. 2b. For Bi < Bq this can lead to interesting 
potentials with sharp repulsive structures, see the inset 
of Fig. 2a. 




FIG. 4: Tunneling vs. decoherence rates. Tunneling 
rate t (left axis) for a square lattice at sufficiently large A 
such that the lattice potential can be approximated by Mat = 
Vb(sin^ kx + sin^ ky), see text. This is plotted for Rubidium 
(solid lines) and Lithium (dashed line) for A = 100 nm, figure 
a), and for a more challenging value A = 20 nm, figure b). 
On the right axis the decoherence rate F = Fsf + Tq^i due to 
thermal jiggling of vortices, see text, is plotted for Rubidium 
(solid line) and Lithium (dashed line). We have used zq — 
dm = and chosen Bm such that cjl = 100 cjt (so that one 
can fulfill cjt <^ (^d ^ cjl), which corresponds to fields up to 
- 1 mT (- 10 mT) for A 100 nm (A 20 nm), /x^^ /xb 
(Bohr's magneton), r]/A = 10~^ N/ms, and a temperature 
of T = 4 K. For a = 1.5 A, one has that for A = 100 nm, 
ujt/27v - 9 X 10^ Hz (- 2 X 10^ Hz) and t/27v - 2 x 10^ Hz 



Decoherence due to thermal jiggling of vortices 



(^ 3 X 10^ Hz) for Rubidium (Lithium). Using A 
then cjt/27r - 2 x 10^ Hz (- 5 x 10^ Hz) and t/27v ' 
Hz (- 9 X 10^ Hz). 



: 20 nm, 
5 X 10^ 



Atoms in magnetic traps are subjected to decoher- 
ence due to magnetic field fluctuations at the Larmor 
frequency cjl- These fluctuations yield spin- flips to un- 
trapped atomic states, thereby leading to atom losses, as 
well as to motional heating [38]. In metal surfaces, these 
fluctuations are generated by thermally excited motion 
of electrons (Johnson Noise), as it is well understood in 
the context of atom chips [39]. Superconducting vortex- 
free surfaces have been predicted to dramatically reduce 
Johnson noise by 6-12 orders of magnitude [40-42]. How- 
ever, experiments in superconducting atom chips have 
shown only a moderate improvement since they operate 
in a regime where uncontrolled superconducting vortices 
are present [43-46]. In the SVL proposed here, it is clear 
that a source of magnetic field fluctuations will be given 
by the thermal jiggling of the pinned vortices. In the 
SI we use a standard phenomenological model of vor- 
tex dynamics (see [48] and reference therein) in order to 
estimate the spin-flip rate Fgf and the motional heating 
rate Fq^i [38] induced by the thermal motion of vortices. 



They are given (up to some constant factors) by 



37r^ 



^0 



P 

kBT 



(6) 



k^ 



UJt 



Here xq = ^/h/{2mg,uJt)^ k^ = ^Q/(2/ioAa^) (provided 
27rA a) is the spring constant given by the repulsive 
force with the lattice, see [47] and SI, and Ud = kp/r] 
with T] being the vortex viscosity coefficient, is the so- 
called depinning frequency, which marks the crossover 
between elastic motion, dominant at lower frequencies, 
and purely dissipative motion, arising at higher frequen- 
cies [48]. Using typical numbers for the vortex viscosity 
77/A = 10~^ Kg/(s m) [48], these rates are remarkably 
small compared to the tunneling rate in Hubbard Hamil- 
tonians, see Fig. 4. In the SI, we also estimate the damp- 
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ing rate in the harmonic motion of the atom due to the 
magnetic couphng to the overdamped vortex, which turns 
out to be neghgible. 

In the actual experiment there might be other sources 
of decoherence that we are not considering here. Prac- 
tical considerations might also be very relevant. For in- 
stance, as analyzed in the SI, the position of the vortices 
has to be very accurate, with an error less than 1 — 2 % in 
the distance between them. Otherwise, the trap depths 
and thereby the tunneling rates will fluctuate through 
the lattice. This can constitute a serious challenge in the 
nanofabrication of regular antidot lattices. In this re- 
spect, using triangular lattices spontaneously formed in 
a film without artificial pinning might be advantageous. 
The randomness in size and shapes of the antidots might 
also lead to imperfections, nevertheless, note that the 
flux in each vortex is given by the constant of nature ^q. 
Note that time dependent fields might induce dissipation 
in the vortices. As discussed before, we apply a time 
dependent field Bm(^) to avoid Majorana losses. This 
field is parallel to the film, and in the ideal case, would 
not interact with the vortices. In any case, to reduce 
dissipation it will be convenient is to use ujm <^ ^d- 

Magnetic local addressing 

Finally, let us discuss the possibility to perform mag- 
netic local addressing of the atoms in the lattice. We will 
use the fact thin films of d ^ X are effectively transparent 
to external magnetic fields, as explained below. Conse- 
quently, we propose to place a magnetic tip close to the 
bottom side of the film to locally interact with the atoms 
trapped above, see Fig. 5a. In the SI, we obtain the an- 
alytical expression of the magnetic field above the film 
for a given thickness d and London penetration depth A. 
The ratio between the magnetic field dit z = zq -\- d in 
the presence of a superconducting film with London pen- 
etration depth A with the corresponding one in case of 
not having the film q depends on the dimensionless pa- 
rameters X/d and L = (zo+ad)/G^, see SI. In Fig. 5b, this 
ratio is plotted as a function of X/d for different L. Note 
that even for A = (i, 5| ~ 5| q/2 for L = 2. As discussed 
in [50] the minimum distance for which the dipole will 
create a vortex is given by ai = A^ /io^d/[*^oA ln(A/^)], 
this equation is valid for ai > A. With typical numbers 
(A = 100 nm and ^ = 10 nm), ai ^ 2 A for a dipole 
of magnetic moment rrid ~ 10^/iB- Using the maximum 
magnetic moment and X = d = A = zq = ad/2 = 100 nm 
(note that L = 3) this leads to a coupling to the atom 
of Qd - I-IbBI/H - 0.4/iB5d,o/^ - 27r X 10^ Hz. Since 
the neighbor atoms are farther away from the dipole, the 
coupling is reduced at least by a factor of two, which 
permits to perform local addressing by adding different 
phases into the internal state of the atoms. In order to 
measure the collective state of the atoms one can release 
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FIG. 5: Magnetic local addressing, a) An atom trapped 
at zq above a film of thickness d is magnetically addressed by 
approaching a magnetic tip at a distance aa from below the 
film, b) The ratio between the field in the presence of the film 
over the same field without it B^/B^^q is plotted as a function 
of X/d for L = (^0 + ad)/d = 1, 2, 4. Note that we recover the 
expected limits: when X ^ d the superconductor completely 
shields the region above the film from the field of the dipole; 
when X^ d the superconductor becomes almost transparent 
to the dipole field. 



them from the trap and perform time-of-fiight measure- 
ments similarly to what is done in optical lattices [1]. 
Recall that the magnetic lattice can be switched off by 
decreasing the bias field Bi. 

Outlook 

To conclude, we have shown that the control and 
manipulation of superconducting vortices in thin films 
can be used to design a magnetic nanolattice for ultra- 
cold atoms. This is based on trapping neutral atoms 
near a surface in a superconducting state, whose macro- 
scopic coherence leads to promising scalings regarding 
surface decoherence. In the quest of bringing quantum 
simulation into a non-classical regime, for instance, by 
simulating antiferromagnetism in Hubbard models and 
long-time quantum many-body dynamics, the SVL pro- 
posed here may provide a very useful path by boosting 
the energy-scales of the simulator by scaling down the 
lattice. The interplay between superconductivity and 
atomic physics, as proposed here, might pave the way 
towards all-magnetic schemes for quantum computation 
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and simulation with neutral atoms as well as to side ap- 
plications such as using cold atoms to probe important 
properties of high-Tc superconductivity. 
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